Effective spinless fermions in the strong coupling Kondo model 
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Starting from the two-orbital Kondo-lattice model with classical t2 g spins, an effective spinless 
fermion model is derived for strong Hund coupling Jh with a projection technique. The model 
is studied by Monte Carlo simulations and analytically using a uniform hopping approximation. 
The results for the spinless fermion model are in remarkable agreement with those of the original 
Kondo-lattice model, independent of the carrier concentration, and even for moderate Hund coupling 
Jh- Phase separation, the phase diagram in uniform hopping approximation, as well as spectral 
properties including the formation of a pseudo-gap are discussed for both the Kondo-lattice and the 
effective spinless fermion model in one and three dimensions. 
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I. INTRODUCTION 



The study of manganese oxides such as Lai-xSr^MnOa 
(LSMO) and Lai^Ca^MnOa (LCMO) has attracted 
considerable attention since the discovery of colossal 
magnetoresistance in these compounds JjJ || . These ma- 
terials crystallize in the perovskite-type lattice structure 
where the crystal field partially lifts the degeneracy of the 
manganese d-states. The energetically favorable three- 
fold degenerate t 2g levels are populated with localized 
electrons, which according to Hund's rule form localized 
S = 3/2 spins. The electronic configuration of the Mn 3+ 
ions is i|g e g with one electron in the e g orbital, which is 
missing in the Mn 4+ ions. The e g electrons can move be- 
tween neighboring Mn ions mediated by bridging 2 ~ 2p 
orbitals. The interplay of electronic, spin, and orbital de- 
grees of freedom along with the mutual interactions, such 
as the strong Hund coupling Jh of the itinerant electron 
to localized t2 g spins, Coulomb correlations, and electron- 
phonon coupling leads to a rich phase diagram includ- 
ing antifcrromagnetic insulating, ferromagnetic metallic, 
and charge ordered domains || . Charge carriers moving 
in the spin and orbital background show interesting dy- 
namical features [Q, |) . The electronic degrees of freedom 
are generally treated by a Kondo-lattice model, which in 
the strong Hund coupling limit is commonly referred to 
as double-exchange (DE) model, a term first coined by 
Zener§. 

Monte Carlo (MC) simulations have contributed signif- 
icantly to our understanding of the manganites. Intense 
MC simulations for the DE model have been performed 
by Dagotto et al and Furukawa || in the space of 
the classically treated t 2g spins. Static and dynamical 
observables of the Kondo model have been determined 
||. These MC simulations gave first theoretical indica- 
tions of phase separation (PS) [[to) in manganite models. 
Preliminary studies have been performed to analyze the 



importance of nearest neighbor Coulomb repulsion in the 
two-orbital DE model [0] as well as the importance of 



classical phonons |12 
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Many publications are based on the Jh = oo limit. 
Here we propose an effective spinless fermion model for 
the strong coupling limit of the Kondo-lattice Hamilto- 
nian, which is equally simple as the Jh — oo limit but 
which still contains the crucial physical ingredients of fi- 
nite Jh- The dynamic variables are e g electrons with 
spins parallel to the t2 g spins at the respective sites. The 
influence of antiparallcl spins is accounted for by the ef- 
fective Hamiltonian. The derivation of the model is based 
on a projection technique, analogous to the derivation of 
the t J model from the Hubbard model. The role of the 
Hubbard U is played by Jh which couples to the classical 
t 2g spins. In contrast to the t J model, the high-energy 
subspace is thus controlled by classical variables and con- 
sequently the resulting model is much simpler than the 
t J model. For a given t 2g spin configuration, the result- 
ing Hamiltonian is a one-particle operator. Its electronic 
trace can be evaluated analytically, once the one-particle 
energies are known, leading to an effective action for the 
t 2g spins, which can be simulated by Monte Carlo tech- 
niques. 

The obvious advantage of this approach is the reduc- 
tion of the dimension of the Hilbert space. This can be 
exploited in MC simulations by going to larger systems 
and/or additional degrees of freedom. 

For a large range of parameters, the effective spinless 
fermion model is found to yield very satisfactory results 
and to perform much better than the rough Jh — > oo ap- 
proximation. We compare spin- and charge-correlations 
as well as quasi-particle spectra of the projection ap- 
proach with the full two-spin model and with the Jh — > 
so limit. 

The effective model is treated without approximations 
by Monte Carlo simulations as well as by a uniform hop- 
ping approximation (UHA) capturing the essential in- 
fluence of the ti g spins on the e g electrons. The UHA 
computation can be performed analytically, particularly 
in the thermodynamic limit. Most of the UHA results are 
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found to be in striking agreement with MC results. We 
find two phase transitions as a function of the chemical 
potential, one close to the empty band and the second 
close to a completely filled band. At each phase transi- 
tion we observe PS, as reported for the upper transition 
in fl(i|| . For a ID-chain we derive an analytical expres- 
sion for the two critical chemical potentials at which (PS) 
occurs. 

For the 3D Kondo-lattice model canonical UHA re- 
sults yield a phase diagram which displays various types 
of antifcrromagnetic (AF) order including spin-canting, 
as well as ferromagnetism (FM). Our finite Jh results for 
3D are in close agreement with those derived in the limit 
of infinite Hund coupling [fL3| . In the grand canonical 
ensemble we find, however, that only the 3D antiferro- 
magnetic and the 3D ferromagnetic order prevail. The 
transition is again accompanied by PS. 

The paper is organized as follows. In Sec. [n]the Kondo- 
lattice model is introduced. By applying a projection 



Finally, in Sec. VI we summarize the key conclusions. 



technique in Sec. Ill this model is mapped onto the ef- 
fective spinless fermion model. In Sec. IV we present 



the phase diagrams and phase separation boundaries in 
one and three dimensions within a uniform hopping ap- 
proximation. Results of Monte Carlo simulations for the 
original and the projected model are discussed in Sec. [v|. 



II. MODEL HAMILTONIAN 

In this paper, we will concentrate on purely electronic 
{tig-i^g) properties, leaving phonon degrees of freedom 
for further studies. As proposed by Dagotto et al. Q| 
and Furukawa || , the ti g spins S, are treated classically, 
which is equivalent to the limit S ^ oo. The spin de- 
grees of freedom are therefore replaced by unit vectors S; , 
parameterized by polar and azimuthal angles 9i and <pi, 
respectively, that represent the direction of the ti g spin 
at lattice site x 4 . The magnitude of the spin is absorbed 
into the exchange couplings. It is expedient to use the 
individual t2 g spin directions as local quantization axes 
for the spin of the itinerant e g electrons at the respec- 
tive sites. This representation is particularly useful for 
the Jh —* oo limit, but also for the projection technique, 
which takes spin-flip processes for finite Hund coupling 
into account. 

It is commonly believed that the electronic degrees 
of freedom are well described by a multi-orbital Kondo- 
lattice model 



H = 
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It consists of a kinetic term with modified transfer inte- 
grals t°£jp, where are site-, a(f3) orbital-, and cr(a') 
spin indices. The number of lattice sites will be denoted 
by L and the number of orbitals per site by M . The oper- 
ators cl aa (c iaa ) create (annihilate) e g -electrons at site Xi 
in the orbital a with spin parallel (a =t) or anti-parallel 
(cr =1) to the local ti g spin orientation S^. The next term 
describes the Hund coupling with exchange integral Jh- 
As usual, hi aa is the spin-resolved occupation number op- 
erator. Usually, the Kondo-lattice Hamiltonian contains 
an additional term proportional to the electron number 
N e , H c = —JhN c , which has been omitted in Eq. ([!]), as 
it merely results in a trivial shift of the chemical potential 
fj,^fi-J H . 

The modified hopping integrals depend upon the 

geometry of the e g -orbitals and the relative orientation 
of the t2 g spins: 



ia,jP 



ia,j/3 



The first factor on the RHS is given by the hopping am- 
plitudes Uajp which read 
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as matrices in the orbital indices a, (3 = 1(2), correspond- 
ing to the x 2 — y 2 (3z 2 — r 2 ) orbitals (see e.g. [Q). The 
overall hopping strength is t, which will be used as unit 
of energy, by setting t = 1. The relative orientation of 
the ti g spins at site i and j enters via 
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(3) 



with the abbreviations Cj — cos(9j/2) and Sj = sin(8j/2) 
and the restriction < 6j < ir. The modified hopping 

part of the Hamiltonian is still hermitian, since = 

Finally, Eq. (nf) contains a super-exchange term. The 
value of the exchange coupling is J' ~ 0.02 0], account- 
ing for the weak antiferromagnetic coupling of the ti g 
electrons. Here we will approximate the local ti g spins 
classically. For strong Hund coupling Jh ^> 1 the elec- 
tronic density of states (dos) consists essentially of two 
sub-bands, a lower and an upper 'Kondo-band', split by 
approximately 2Jh- In the lower band the itinerant e g 
electrons move such that their spins are predominantly 
parallel to the ti g spins, while the opposite is true for the 
upper bandp4[. Throughout this paper, the electronic 
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density n (number of electrons per site) will be restricted 
to < n < 1, i.e. predominantly the lower Kondo-band 
is involved. 

III. PROJECTION TECHNIQUE 

The separation of energy scales M| is a well known 
strategy to simplify quantum-mechanical many-body sys- 
tems. In the case of manganites, the Hund coupling Jh 



is known to be much greater than the other parameters t 
and J'. Consequently, the hopping to antiparallcl e g — t 2g 
configurations can be treated in second order perturba- 
tion theory [TtJ by a projection approach. On the 
low energy scale, the dynamical variables are e g electrons 
with spin parallel to the local ti g spins. The virtual ex- 
citations (ia t) — > (j/3 I) — > (if a' j), which are medi- 
ated by the hopping matrix, lead to an effective spinless 
fermion Hamiltonian 
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The effective Hamiltonian contains the kinetic energy of 
e g electrons with spin parallel to the ti g spins (first term). 
The kinetic energy is optimized by aligning all tig spins 
which is the usual ferromagnetic double exchange effect. 
The second term describes an additional hybridization 
and favors antiferromagnetic t2 g spins leading to an effec- 
tive antiferromagnetic interaction J e ff which is generally 
stronger than J'. The "three-site" hopping processes of 
the third term are of minor influence. We will see that 
this term is in general negligible. On the other hand, 
its inclusion does not really increase the numerical effort. 
Eq. (||) is valid for arbitrary hopping matrices U a ^p. In 
subsequent sections, however, the discussion will be re- 
stricted to nearest-neighbor hopping only. 

The Hamiltonian Eq. (EJ) constitutes a spinless fermion 
model, similar to the one obtained in the Jh — > oo-limit, 
which can be treated numerically along the lines pro- 
posed by Dagotto and coworkers and Furukawa ||. 
Finite Jh values can thus be treated with the same nu- 
merical effort as the case Jh = oo. In the MC simulations 
the weight for a ti g spin configuration is determined by 
the grand canonical trace over the fermionic degrees of 
freedom in the one-electron potential created by the ti g 
spins. 

The obvious advantage of Eq. (Q) as compared to 
Eq. (|l|) is the reduced Hilbert space. 



IV. UNIFORM HOPPING APPROXIMATION 

Before discussing approximation-free MC results for 
the effective spinless fermion model, we will investigate 
the main features of the Hamiltonian (Q) by a uniform 
hopping approximation proposed by van-den-Brink and 
Khomskiilla]. To this end, we introduce two different 
mean angles between neighboring ti g spins, one in z- 
direction (6 Z ) and one in the xy-pl&ne (0 xy ). It should 
be stressed that Z and 6 xy are relative angles between 



adjacent spins, with values between and tt, and are 
not to be confused with the polar angles 0j. We as- 
sume that these angles are the same between all neigh- 
bor spins, i.e. Sj • Si± z = cos9 z for all lattice sites i 
and Si ■ Si±x — Si ■ Si± g = cos9 xy . The allowed spin 
configurations include, among others, ferro- and antifcr- 
romagnetism as well as spin canted states. The impact 
of the t2 g spins on the hopping amplitudes simplifies to 

nr = cos(|) , = sin(|). 

for hopping processes along the z-direction and similarly 
for u^y for electron motion in the xy-plane. The hop- 
ping matrix is now translationally invariant. The in- 
ner product of the tig spins entering the super-exchange 
reads 

Sj ■ S l+i = cos 8 Z =2u 2 z -l 
and similarly for neighboring pairs in the xy-planc. 



A. Phase Separation in ID systems 

First we consider the simplest case, namely a ID chain 
in z-direction and we ignore the additional three-site hop- 
ping (third term in Eq. (El)). At the end of this section 
we will show that it can indeed be neglected. Due to the 
symmetry of the hopping elements the x 2 — y 2 -orbitals 
form an irrelevant dispersionless band, which will be ig- 
nored in the sequel. The influence of the average spin 
orientation is captured in the uniform hopping ampli- 
tude u. Assuming periodic boundary conditions, the 
Hamiltonian simplifies to 

H = _u E c i°i _ ~~J~ E 4^ + J'L(2u 2 - 1) , (5) 
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where we have dropped orbital indices. The virtual hop- 
ping processes couple merely to the density and the dis- 
persion of the spinless fermions is given by the shifted 
tight-binding band structure 



e fc = -2ucos(fc) - (1 - u 2 )/J H 



(0) 



The band width is 4u. It vanishes accordingly for AF or- 
der and reaches a maximum for FM order. In Fig. |l| the 
resulting band-filling is schematically depicted for zero 
temperature as function of chemical potential and hop- 
ping amplitude. The condition for an empty/ filled band 
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FIG. 1: Filling of the tight-binding band depending upon 
chemical potential and hopping amplitude. Condition for par- 
tial filling: u > u*(fi). The minimum free energy solution in 
UHA (dashed line) exhibits a jump (phase separation) from 
AF to FM order at /xf and back to AF order at 

depends on the 'effective chemical potential' 

f i :=( fI +J-(l-u 2 ))/u. (7) 

According to /i < min^e^), the band is empty if /2 < —2. 
For the completely filled band the condition reads fi > 2. 
Partial filling is possible for intermediate values of the 
chemical potential (—2 < \i < 2) if the hopping ampli- 
tude exceeds a threshold it*(/x). The logarithm of the 
grand canonical partition function reads 



hxY 



] )-f3J'L(2u 2 - I) 



In the thermodynamic limit (L 
free energy per lattice site is 



oo) and for T — the 



/ = uE{p) -ujl N{p) + J'(2u 2 - 1) 



(8) 



with E(x) being the mean kinetic energy and N(x) the 
mean particle number of a tight-binding band with dis- 
persion — 2cos(fc). For \x\ < 2 these quantities are 



E(x) = 



V4T 



, . 1 arcsin(x/2) 
N(x) = - + 

Z 7T 



The kinetic energy E(x) is zero for the empty band 
(x < —2), as well as for the completely filled band 
(x > 2). The mean particle number N(x) is zero if the 
band is empty and unity if it is full. Fig. ^| shows the free 




FIG. 2: Free energy as a function of chemical potential and 
hopping amplitude u for T — 0, Jh = 6 and J' = 0.02. The 
solid lines are for ^1 and ii2, respectively. For the sake of 
clarity, f(u = 0, fi) has been subtracted. 



energy as a function of chemical potential and hopping 
amplitude u. We find local minima at u = (AF order) 
and u = 1 (FM order). The kinetic terms decrease with 
increasing u, favoring FM order, while the ti g spin ener- 
gies increase with increasing u, favoring AF order. The 
global minimum switches from AF to FM at the critical 
chemical potential /i — [i\ and back to AF at /i = 
The values for the critical chemical potential follow from 
the condition /|„ =0 = /| u =i, yielding 

(// + 1/J H ) 7V (/i c ) + 2 J' - v c N(p c ) - E(p c ) , (9) 

where No(fi c ) denotes the mean particle number for 
u = 0, i.e. for perfect AF order. In this case, the tight- 
binding band is dispersionless and Nq(^ c ) is either zero 
or one, depending upon the actual value of the chemical 
potential. 

For the standard parameter set Jh = 6 and J' = 0.02 
the numerical values for [i c are [i\ = —1.6730 and = 
1.03431. 

The UHA solution corresponds to the global minimum 
of the free energy. Its location is depicted in Fig. [I] and 
the corresponding densities are shown as lines in Fig. |^. 
For large negative chemical potential the system is an- 
tiferromagnetic and the e g band is empty. At /i^, AFM 
domains with zero electron density coexist with FM do- 
mains with finite density n\. Increasing \i leads to fer- 
romagnetism, and the filling increases gradually with /x, 
following the tight binding formula n = arccos(— fi/2)/n. 

At fi2, FM domains with density n-i > n\ coexist with 
AFM domains of density one. Finally, for /i > /zJj, the 
system jumps back to antifcrromagnetism, now at density 
one. We thus see that the system exhibits phase separa- 
tion. It should be pointed out that PS is suppressed if 
nearest-neighbor Coulomb repulsion among e g electrons 
is included into the model p8|. 

Let us now discuss the values of n\ 2 and the size of 



5 




-12 -10 -8 -6 -4 -2 



H-J H 

FIG. 3: Electron density versus chemical potential for Jh = 
4,6, and 10 (right to left), and J' = 0.02. The lines corre- 
spond to UHA. MC results at /3 = 50, L = 20 for the spinless 
fermion model H v (circles) are compared with those for the 
DE model H (crosses). Error bars of the MC data are smaller 
than the symbols. 
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FIG. 5: Discontinuity in n as a function of the effective ex- 
change coupling Jeff. 

sion of Eq. (0) with respect to u c about u c = — 2 yields 
An = (3J' /tt) 1 / 3 . At Jeff = the slope of the curve di- 
verges, implying that already an infinitesimal J e ff leads 
to phase separation in this regime. 

For realistic parameter values (Jh = 6 and J' = 0.02) 
we find Ani 0.18 and An2 ~ 0.3, respectively. The 
second value is mainly driven by the virtual hopping, 
which increases the tendency towards PS. 



B. Impact of "three-site-terms" 



0.5^ 1 1 1 

0.05 0.1 0.15 

FIG. 4: Dependence of the critical chemical potential 1/^2 
on the effective exchange coupling J c g. 

the discontinuity in n. According to Eq. (^), the first 
critical chemical potential n\ = p c (J'), corresponding 
to Nq(/i) — 0, is independent of the actual value of 
Jh- Here the effective antiferromagnetic interaction J e s 
purely stems from the superexchange coupling of the tig 
spins, J e ff = J'. This does, however, not mean that the 
Hund coupling is irrelevant for this phase transition. On 
the contrary, the phase transition is driven by the FM 
tendency introduced by the Hund coupling. The inde- 
pendence of Jh in our results means that there are no 
second order corrections to the Jh — oo limit. The de- 
pendence of /i c on J' is depicted in Fig. |^. 

Now we turn to the second phase transition corre- 
sponding to Ao(/i c ) = 1. This transition is controlled 
by the stronger effective exchange coupling 

J eff = J ' + > ( 10 ) 

as can be seen from Eq. (^|), or directly from Eq. (||) 
at N = 1. Due to the particle-hole relation N(—[j,) + 
N(fi) — 1 the second critical chemical potential is given 
by Ml = ^^ c {Jcs), depicted in Fig. |[ In the limit 
Jh — * oo we have \i\ = — The density discontinuity 
is An := N(p!{) = An(J') at n\ and An := 1 - N(nQ = 
An(J G ff) at It is depicted in Fig. || Series expan- 



Here we consider the impact of the additional hopping 
in Eq. ([|), which in the ID case with one orbital results 
in a next-nearest neighbor hopping 

1 — ii? 

i 

Combined with the terms of Eq. (JsJ) , the resulting single- 
particle dispersion reads 

1 -u 2 l-u 2 

efc = —2m cos k — — — — cos(2fc) . 

Jh Jh 

In the limit u — > 1 we recover the original tight-binding 
band. The density of states has additional van Hove sin- 
gularities. Contrary to the dispersion of Eq. (||), the band 
width remains finite in the limit u — > 0, due to H*. 

Next we derive the conditions for phase separation. In 
the limit T — > 0, the free energy is given by 

/= T dcp„(e)(e-M) - J'(2m 2 -1) , 

J — oo 

where p u (e) denotes the density of states corresponding 
to £fc . Numerical evaluation shows that the minima of / 
are still at u = and u = 1. The condition for phase 
separation is therefore still /| u =o = f\u=i- In principle, 
due to the finite width of the band at u — intermediate 
particle numbers iV (u c ) are possible. A detailed calcu- 
lation shows, however, that for realistic parameters, only 
No(/j, c ) = or No(ii c ) = 1 can meet the phase separa- 
tion criterion. For Nq(/i c ) = and Nq(ii c ) = 1, however, 
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no hopping is possible and consequently the additional 

hopping term vanishes. Therefore, the criterion for PS is 

the same as in Eq. (^) and the three-site hopping has no 0.02 

influence on the critical chemical potentials /i c at which 

phase separation occurs. 0.04 

In general, the modification of the bandwidth due to 
next-nearest neighbor hopping is small. It has almost q qq 

negligible impact on the results. On the other hand it 
poses little extra-effort to include it in MC simulation. . no 

O.Oo 



C. Phase Diagram in 3D systems 



0.1. 



0.1 



0.2 



0.3 



0.4 



0.5 



We now turn to the 3D case. In uniform hopping ap- 
proximation, the super-exchange reads 

H se = J'L 3 (2(2u 2 xy -l) + (2u 2 z -l)) 

where L denotes the linear dimensions of the sys- 
tem. Upon substituting uniform hopping amplitudes into 
Eq. (0) , the fermionic part of the Hamiltonian can easily 
be diagonalized. The one-particle energies are given by 
the eigenvalues of a 2 x 2 matrix with matrix elements 

£ n( fc ) - u xy (cosk x + cosfcy) - (1 - u 2 xy ))/J H ^ 

ei2(fc) =— u xy (cosk x - cosfc y ) = e 2 i(fc) 

^22 (k) = — 2u z cos k z — u xy ( cos k x + cos k y )/2 
-(l-u 2 z + (l-ul y )/2)/J H 

where the subscript 1 (2) refers to x 2 — y 2 (3z 2 — r 2 ) or- 
bitals. The symmetry of the e ff -wavefunction has been 
exploited in the above expressions. As a consequence of 
the UHA with two different angles 8 Z and 9 xy , virtual 
hopping processes cannot induce transitions between dif- 
ferent orbitals, and Jh appears only in the diagonal ele- 
ments of the matrix. 

We determine the phase diagram in the canonical en- 
semble at T = with respect to electron density n and 
exchange coupling J' for fixed Hund coupling Jh = 8, 
with Eq. ( |TT| ) evaluated on a 20 3 momentum lattice. For 
each parameter set, the free energy is minimized with 
respect to the hopping amplitudes. The resulting phase 
diagram is depicted in Fig. ^. At very low doping, 3D an- 
tiferromagnetic (G) order dominates, irrespective of the 
value of J', as long as J' > 0. Increasing the electron 
concentration for J' > 0.02, the system favors first a C- 
phase, then an A-phase, and finally ferromagnetism (F) . 
Similar results have been found for the Jh = oo limit p3[ . 
Finite values for Jh have almost negligible influence on 
the phase diagram for densities n < 0.5. The G-phase has 
not been reported in Jl^| since it has not been taken into 
consideration. For small super-exchange of the t2 g spins 
(J' < 0.018) the transition from G to A phase evolves via 
spin canting. By increasing the electron doping, the F- 
phase is reached without canting. The situation is more 
complex for larger values of J'. Fig. ^ shows the optimal 



FIG. 6: Canonical T — UHA phase diagram of the spin- 
less two-orbital Kondo-lattice model with classical t^g spins 
for Jh = 8. Depending upon the electron concentration of 
the canonical ensemble and the super-exchange, G-type, A- 
type, Gtype antiferromagnetic phases or the ferromagnetic 
phase (F) are observed. Meaning of phases: G=(AF,AF); 
A=(FM,AF); C=(AF,FM); F=(FM,FM), where the first en- 
try denotes the order in the xy-plane and the second in 
z-direction. The solid line represents the phase-boundary 
between G- and F-order, respectively, obtained in a grand 
canonical ensemble. 



angles 6 Z and 9 xy as a function of the electron density n 
for J' = 0.025 and J H = 8. For n < 0.039, there is 3D 
AF order (G-type) with an increasing tendency towards 
canting between ti g spins in the xy plane. At n ~ 0.04 
this tendency is strongly reduced while at the same time 
cos(# z ) discontinuously jumps to zero and we gradually 
enter the C-phase by aligning the spins along the chains 
in ^-direction. The ferromagnetic chains are not perfectly 
antiferromagnetically stacked. At about 8% electron den- 
sity we observe a phase transition to the A-phase. At 
n ~ 0.15 an abrupt transition to 3D ferromagnetic order 
occurs. 




FIG. 7: Evolution of the optimal angles 6 Z 
of the electron density n at J' = 0.025, Jh = 8. Further 
details like in Fig. ^. 



Besides the analysis for the canonical ensemble, com- 
putations have been performed for the grand canonical 
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ensemble as well. The results are significantly different. 
In the grand canonical ensemble only the G and the F 
phase remain. The white solid line in Fig. || represents 
the phase boundary between the two phases. For fixed 
J', the behavior is similar to that of the ID system, de- 
picted in Fig. |[ Below a critical chemical potential /xf 
the electron density is zero and the ti g spins have AF 
order. At fi — zero density and a finite density ni, 
given by the solid line in Fig. |[ coexist. Concurrently 
with phase separation, AF- and FM order coexist. Above 
n.\ the density increases monotonically. The second tran- 
sition at /iJ; is not shown in Fig. ^ as it occurs close to 
n = 1. Therefore, the grand canonical UHA result does 
not exhibit the additional magnetic phases (A and C), 
which are observed in experiment. The relevant densi- 
ties are never stabilized. 



V. MC SIMULATIONS FOR ID SYSTEMS 

In this section we compare MC results, obtained for 
the original double exchange (DE) model (|l|), with those 
for the effective spinless fermion model (Eq. (|J)), where 
the additional ("three-site") hopping term has been ne- 
glected. We use the grand canonical Monte Carlo method 
introduced in JhJ with open boundary conditions (obc). 

We restrict the discussion to ID systems. There is no 
reason to believe that the performance of the approxi- 
mation of effective spinless fermions is different in higher 
dimensions. Furthermore, the approximation has little 
influence on the orbital degrees of freedom and we restrict 
the analysis to the one-orbital model. In all simulations, 
the super-exchange coupling of the ti g spins is J' = 0.02. 

Fig. H shows the dependence of the electron density on 
the chemical potential. The system parameters are L = 
20, (5 = 50, and Jh is varied between J H — A and Jh = 
10. All results for the two models are in almost perfect 
agreement. The 'largest' difference can be observed at 
for Jh = 4. The lines in Fig. || represent the results of 
the uniform hopping approximation, which are strikingly 
close to the MC data points. The discontinuities are more 
pronounced in UHA than in the MC data, which can 
partially be attributed to the fact that the UHA results 
are for L = 00 and T = 0. The treatment of finite 
temperatures in the UHA requires the determination of 
the number of ti g configurations at a given u and is the 
subject of current investigations pq . 

The structure factor of the ti g spins has been calcu- 
lated for various densities n in the grand canonical en- 
semble by adjusting the chemical potential. The results 
are illustrated in Fig. ||. Again the data for the two mod- 
els, Kondo- and effective spinless fermion model, are in 
perfect agreement within the error bars. Corroborating 
the UHA results of the previous section, the filled band 
(?i = 1) has a peak at k = 7r, corresponding to AF order. 
For decreasing density the ferromagnetic peak increases 
up to n = 1/2 and then it decreases again. The inset 
shows the nearest neighbor and next nearest neighbor 



spin correlation function versus density. We observe that 
both models yield the same magnetic behavior: AF or- 
der at low and high density and a ferromagnetic phase 
at intermediate fillings. The pronounced peak at k = tt 
for n — 1 results from (virtual) spin-flip processes, driven 
by the relatively strong exchange coupling J c ff = 0.103, 
Eq. (10). Contrariwise, the AF structure near n = is 
much less pronounced as it is merely driven by the weak 
super-exchange coupling J' = 0.02. Fig. KJ shows the 




FIG. 8: Structure factor for the t2 a spins at various electron 
densities with Jh = 6, /3 = 50, and L = 20, for the spinless 
fermion model H p (circles) , and for the DE model H (crosses) . 
When not shown, error bars are smaller than the symbols. 
The inset shows the nearest neighbor (solid line) and next 
nearest neighbor (dashed line) t2 g spin correlation. 

structure factor for the ti g spins at n = 1 for different 
values of Jh- The AF peak decreases with increasing 
Jh and degenerates to a broad structure in the limit 
Jh — ► 00. Obviously, the inclusion of the second or- 
der term to the effective spinless fermion model, which is 
missing in the commonly used Jh = 00 limit, and which 
provides the strong exchange coupling J G ff, is crucial for 
the correct description of the AF order at high electron 
density. The inset shows the spin structure for n « 3/4, 
at which the system exhibits ferromagnetic order. In this 
case, the FM correlations increase with increasing Jh- 

As a further test for the spinless fermion model, the 
electronic contribution to the total energy is shown in 
Fig. [h]. The energy for the Kondo-lattice model (first 
two terms in Eq. ([!])) is compared with the correspond- 
ing contributions in the effective spinless fermion model 
(Eq. (^)). For all fillings the results are in very good 
agreement. 

Quantitatively, the largest differences are found for 
n = 1. For this density the dependence of the energy on 
Jh is studied in Fig. [Fy. A detailed comparison reveals 
that the effective model describes the electronic energy 
extremely well, even in the moderate coupling regime. 
The ubiquitous Jh = 00 approximation, on the other 
hand, yields zero electronic energy. For the parameters 
underlying Fig. [H], the ti g spins are antiferromagnetically 
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FIG. 9: Spin structure factor at n = 1 (inset: n ~ 0.75) for 
different values of Jh- Same symbols and parameters as in 
Fig. H. In the limit Jh — > oo (dashed line) the intensity of 
the AF peak is considerably reduced. 
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FIG. 10: Electronic contribution to the total energy versus 
electron number for = 50, L = 20, and Jh = 6. Results for 
the DE model H (crosses) are compared with those for the 
spinless fermion model H p (circles). 



that the additional (three-site) hopping term in Eq. (0) 
has no influence on this result as the band is entirely 
filled. 

Next we study the spectral function A^iui) in the grand 
canonical ensemble for various mean electron densities 
covering the regimes for AF and FM order, as well as 
phase separation. In all cases the system geometry is a 
20-site chain with open boundary conditions at inverse 
temperature j3 = 50, and exchange couplings J' — 0.02, 
and Jh = 6. We start out with the spectra in Fig. [l^ for 
strong FM order at a mean particle density of n = 1/2. 
According to the inset of Fig. || the spin-correlations are 
0.82 and 0.67 for nearest and next-nearest neighbors, re- 
spectively. The spectral function, depicted in Fig. [1^, 
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FIG. 12: Spectral function obtained by a grand canonical 
simulation for n = 1/2 with Jh — 6, L — 20 and j3 = 50. An 
intrinsic linewidth 7 = 0.1 has been added. The results for 
the Kondo model and the effective spinless fermion model are 
indistinguishable . 
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FIG. 11: Electronic energies versus Hund coupling Jh for 
a completely filled 20-site chain at T = 0. Kinetic energy 
(dotted line) and total electronic energy (dashed lines) for 
the Kondo-lattice model. The solid line represents the total 
electronic energy for the effective spinless fermion model. 



ordered and the lower Kondo-band, or rather the single 
band of the spinless fermion model, is completely filled. 
Nevertheless, the kinetic electronic energy is finite due 
to (virtual) spin-flip processes. It should be pointed out 



resembles closely that of a tight-binding model, valid for 
perfect FM order. The band-width is slightly reduced, 
and for k- values close to the Fermi momentum, A^iu) 
exhibits some minor shoulders. The width (HWHM) of 
the peaks agrees with 7, the value by which the finite- 
size delta peaks have been broadened. The inset displays 
the density of states (dos), which agrees with the tight- 
binding density of states for open boundary conditions. 
Next we increase the mean density to n = 0.75, corre- 
sponding to a chemical potential close to The spin 
order is still predominantly ferromagnetic. The results in 
Fig. [l^ show that both models yield very similar results, 
namely a tight-binding type of quasi-particle band. The 
spectral peaks are, however, significantly broader than 
the mock width 7 and upon approaching the Fermi mo- 
mentum the width increases asymmetrically towards the 
Fermi level. The origin of the broad structures are the 
random deviations of the t^g spins from perfect FM or- 
der. The resulting density of states has piled up spectral 
weight in the center and reveals a precursor of a pseudo- 
gap at the Fermi energy. Interestingly, the dos is almost 
center-symmetric and a mirror image of the 'pseudo-gap' 
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FIG. 13: Same as Fig. |l2|. The indicated particle numbers 
correspond to n = 0.75 n — 0.95, and n = 0.05. Results for 
the Kondo model are represented by solid lines and those for 
the effective spinless lermion model by dashed lines. 



occurs also at the lower band edge. 

The next panel depicts results for n = 0.95, correspond 
to a chemical potential slightly above ^ C 2, where the 
spin order is predominantly antiferromagnetic. Now the 
pseudo-gap, as discussed by Moreo et al. Jl9| , is clearly 
visible. Additionally, we observe a 'mirror pseudo-gap' 
at the lower band edge. There is still good qualitative 



agreement between the results of the two models. Quan- 
titatively, however, there are deviations in the structures 
below the pseudo-gap. The density of states is still re- 
markably well described by the spinless fermion model. 

In the opposite limit of low carrier concentration (n = 
0.05), the chemical potential is close to y, c \ , where we find 
similar features. Again, together with the coexistence 
of FM and AFM order a pseudo-gap shows up at the 
chemical potential as well as at the upper band edge. 
The pseudo-gap is less pronounced in this case, where the 
antiferromagnetic exchange coupling J' is much smaller 
than J e ff at yu.^. 

In both models a considerable amount of spectral 
weight is transferred from the band edges to the cen- 
ter. Interestingly, the density of states is almost center- 
symmetric. The pseudo-gap is present in the spectral 
function irrespective of the wave vector k. Moreo et al. 
argued that the pseudo-gap is formed due to the pres- 
ence of mixed phases with irregular formations of FM 
domains. In contrast to this interpretation, we find the 
pseudo-gap also in the perfect AF regime with a single 
electron (hole). 

Generally we observe, in agreement with ARPES ex- 
periments |2(J , that the width of peaks increases towards 
the Fermi energy and the spectral intensity decreases 
since spectral weight is transferred to the unoccupied 
part of the spectrum, which is not visible in ARPES. 
Furthermore, the peaks are generally much broader than 
the experimental resolution. 



VI. CONCLUSIONS 

We have developed an effective spinless fermion 
model for the strong coupling multi-orbital Kondo-lattice 
model. The effective model has a reduced Hilbert space 
and is particularly suitable for MC simulations. The nu- 
merical complexity is the same as that of the Jh = oo 
model. The reduced Hilbert space allows one to study 
higher spatial dimensions and/or additional degrees of 
freedom, such as phonons. 

Based on the evaluation of various observables the ef- 
fective spinless fermion model performs strikingly well, 
even for moderate Hund coupling. 

It appears that virtual spin-flip processes included in 
our approach, which are missing in the Jh — oo model, 
are crucial for the antiferromagnetic phase close to half- 
filling (n = 1), where they provide a strong effective ex- 
change coupling J c ff = J' + ■ 

Two phase transitions from AF- to FM-order and vice 
versa are observed, accompanied by phase separation. 
Analytic expressions for the chemical potential at which 
phase separation occurs in a ID chain have been derived 
in uniform hopping approximation (UHA). It has been 
shown, however, that they are in extremely good agree- 
ment with approximation-free MC results. 

The UHA phase diagram of the 3D spinless fermion 
model has been determined. In canonical ensembles, the 
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magnetic phase diagram is in qualitative agreement with 
that obtained in previous studies for the Jjj = oo limit. 
Experimentally observed phases, like G-, A-, C-, and F- 
order are found. On the other hand, grand canonical 
ensemble calculations show that only 3D AF and 3D FM 
order prevail. The transition between the two phases is 
accompanied by phase separation. Densities, required for 
other phases, are not stable in UHA. 

The spectral functions show a remarkable center- 
symmetry. In the AF phase, at low and high electron 
density, a pseudo-gap structure is observed at the chem- 
ical potential and a mirror image at the opposite edge of 



the spectrum. 

In passing, it should be noted that nearest-neighbor 
Coulomb repulsion of the e ff -electrons in a two-orbital 
model can be detrimental for phase sep aration. Detailed 
results will be discussed elsewhere 
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